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The purpose of this article is to find an efficient and optimal method 
of compression by reducing the file size while retaining the information 
for a good quahty processing and to produce credible pathological reports, 
based on the extraction of the information characteristics contained 
in medical images. In this article, we proposed a novel medical image 
compression that combines geometric active contour model and quincunx 
wavelet transform. In this method it is necessary to locahze the region 
of interest, where we tried to locahze aU the part that contain the pathological, 
using the level set for an optimal reduction, then we use the quincunx wavelet 
coupled with the set partitioning in hierarchical trees (SPIRT) algorithm. 
After testing several algorithms we noticed that the proposed method gives 
satisfactory results. The comparison of the experimental results is based on 
parameters of evaluation. 
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1. INTRODUCTION 

Nowdayas, with the technological development, the various apphcations such as medical imaging 
require a big space of storage of information or for the transmission, to overcome these constraints, wavelet 
compression is an inevitable solution. Given the importance of medical images and the many pieces 
of information they contain, the accuracy of the most important regions and objects in the image is necessary, 
because it assists an accurate and efficient analysis with credible pathological report with less error. 
Renee the use of wavelet in a separable way does not allow preserving the outlines and the region of interest, 
for this purpose the image segmentation is an independent step. The deformable models allows separating 
image into several regions related with the same properties, and concentrate on specific regions, 
they are defined as a curve or a contour placed in image and deforms to reach an optimal position. 
Recently various methods of medical image segmentation have been explored and presented in 
the hterature [1-3] such as level set [4]. 

The level set (LS) method explained as the operation that aims to gatherthe pixels of image in regions 
that constitutes an image partition according to predefined criteria. The idea behind operation is that an initia 1 
contour specified by the user is moved to the desired objects boundaries by image driven forces. Where we can 
specify two forces, while the deformation process started, the internal forces are designed to conserve the model 
smooth, where the external forces serve to move the model until the object boundary. The aim of our work is 
to use an hybrid method, combines the level set with wavelets, for extracting 
the characteristic information contained in the medical image. Over the years, the wavelets (DWT) have had 
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a huge success in image processing field [5]. The idea of wavelet analysis is to decompose a signal on 
a functional basis of a functional space with weU-defined properties [6]. It is to analyze a signal looking for 
a more compact representation located simultaneously in time and frequency. The wavelet transforms (WT) 
requires less of calculations on account of their inherently multi-resolution nature, and the encoding schemes 
which are suited to applications with important scalability and tolerable degradation 

Eor this paper and including to these two methods and their advantages, we introduce our hybrid 
method by segmenting the desired object, separate the foreground from the background object before applying 
a different wavelet transforms coupled with the SPIHT coding. To properly evaluate 
the compression image quality, we compare the results obtains from various application by the evaluatbn 
parameters. Our paper is organised to give, a general introduction, an explanation of level set (LS) method, 
and quincunx wavelet transform. EinaUy, the results associate to image quality changes 
are presented and discussed. 


2. RELATED WORK 

In recent years, the scientific research has shown an increasing interest in the domaine ROI 
in medicalimage compression, several reasearchers have used many techniques, in order to obtain an optimal 
compression without losing the quality of the information contain the ROI. In this section, recent research 
on ROI compression that have: 

- In 2015, Manpreet and Vikas [7] divided the image into ROI and Non ROI and used fractal lossy 
compression for Non ROI image and Context tree weighting lossless for ROI part of an image. 

- In 2016, Joshi. p andRawat.C [8] segmented the image using maximum entropy based thresholding into 
ROI and Non ROI. After that, they applied lossless compression for ROI using arithmetic coding and lossy 
compression for Non ROI using SPIHT, this method is evaluated in terms OP PSNR, SSIM. 

- In 2017, abdelaziz. E et al [9] introduced a new medical image segmentation approach based on active 
contour improved by ROI extraction. 

- In 2018, C. Priya et al [10] separated the ROI from the Non ROI using fuzzy C means clustering, the Non 
ROI is compressed using CAVLV method and discret wavelet transform compression method is applied 
for ROI. 

- In2019,Nafees A et al [11] compresed the image using DWT based zoning technique in combination with 
DCT for image compression. DWT divides an image into LL, LH, HE and HH frequencies and Zoning 
is further dividing these images into four parts as an input to DCT one after another. The output of DCT on 
each zone is then combined into a compressed bitstream image. 


3. PROPOSED METHOD 
3.1. Level set 

Various techniques and methods are designed and applied for image segmentation, 
are used in the literature [12-13]. It was first introduced by Osher and Sethian [14] who proposed an effective 
implicit representation for evolving curves and surfaces, it allows for automatic change of topology, such as 
merging and breaking. 

The basic idea for image segmentation is to perform the curves as the zero LS of a higher dimensional 
hyper-surface. In the LS form, an interface C is depict implicitly as a LS of a function cp of higher dimension, 
which compute the geometric characteristics and the motion of the front. The interface 
is now represented implicitly as the zero-th LS (or contour). Specific closed curve C, where C=5co, 
with COcQ being an open subset, for the image I(x,y) on the image domain Q. Either the pixel lies on the curve 
itself if (p=0, or it represent the minimum distance between the pixel and the curve. 

According to this equation ^ [15] a closed parametric curve is evolved. 

dt 

Where t is time, E is the rate of evolution andN is the unit normalto the curve. Eigure 1 shows. The curve C 

evolved according to ^ /r/y. Eigure 2 shows level set segmentation steps. 
dt 
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Figure 1. The curve C evolved according to ^ 

dt 


Considering that ^{x^y)> 0 if the point (x, y) inside C, ^(jc,y) < 0 if the point (x, y) is outside C and 
^(jc, y)=0 if the point (x, y) is on C. Thus, the energy functionalF (cl, c2, C) reformulated in terms 
oftheLS function (x, y) as follows [16]: 


/(jv, y) — cl/^ y^^djcdy ( 1 ) 

Cl 

J//(-V, y^ — c:2/^(l — y)y)djcdy 

Cl 

-h //J y))/V^(jv, y)/ dxdy 

Cl 

Where H(Z) and 5d^z) are, respectively, the regularized approximation of Heaviside function H and delta 
function 5 as follows: 


//(Z) = 


1 z > 0 
0 if z-<0 


Siz) = ^Hiz) 
dz 


The LS function (x, y) is updated by the following gradient descent method: 


3(p 

at 




judid^ 






Where cl andc2 are expressed: 


cl(^) = 


J/(x, y^Hsi^ix, y)')dxdy 

Cl 


J Hsijcpi^x^ y^^dxdy 

Cl 


c2(^) = 


J/(x, y)(1 — y^^yixdy 

Cl 

— Hsi(p(^x, yyyydxdy 

Cl 


( 2 ) 


(3) 


(4) 
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(a) 


(b) 


(c) 


Eigure 2. Level set segmentation steps, (a) Curve initialization, (b) Level set segmentation, 

(c) Eoreground extraction 


3.2. Quincunx wavelet 

The separable dyadic analysis require three families of wavelets, which is sometimes regarded 
as a disadvantage, in addition the factor of addition between two successive scales is 4 which may seem high. 
It is possible to solve these two problems, but at the cost of the loss of filter separability and therefoie 
a slightly higher computational complexity. An analysis has been particularly well studied to find a practical 
application, known as "quincunx" [17]. Quincunx decomposition results in fewer subbands than most other 
wavelet decompositions, a feature that may lead to reconstructed images with slightly lower visual quality. 
The method is not used much in practice, but [18] presents results that suggest that quincunx decompositbn 
performs extremely well and may be the best performer in many practical situations. Eigure 3 illustrates this 
type of decomposition. We notice that the dilation factor is not more than 2 between two successive resolutions, 
and that only one wavelet family is necessary [19, 20]. 
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Eigure 3. Quincunx wavelet decomposition 


In this case the Dilatation matrice win be; 


M= 


1 1 

1 -1 
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The Grid transformation (lattice) is done according to Ji -\-l M the following diagram: 


M 


This matrix generates a quincunx lattice in 2D. The column vectors of this matrix form a basis to this 
lattice. The volume of the unit cell associated equals 2. The same lattice Figure 4 is also emanating from 
the matrix below [17]: 

"1 -I 


It is noticed that the dilatation step is ^ ^ on each direction and the geometry of the grid obtained 
justifies the name given to this multiresolution analysis. 



Figure 4. Examples of a lattice quincunx and unit cell, (a) Lattice quincunx, (b) Unit cell 


3.2.1. Quincunx sampling and filter 

First, we recall some basic results on quincunx sampling and perfect reconstruction 
filterbanks [21, 22]. The quincunx sampling lattice is shown in Figure 4. Figure 5(a) shows Quincunx lattice 

and Figure 5(b) shows the corresponding Nyquist area in the frequency domain Letx[/2 ] with n =(ni,n 2 ) G 

denote the discrete signal on the initial grid. The two-dimensional (2-D) z-transform ofx[n] is denoted by: 


x|/i 1^ ? where ^ ^ ^ 


The continuous 2-D Fourier transform is then given by 


(5) 




with 


VF=(vFi,W2) 


( 6 ) 


and finally, the discret 2-D Fourier transform for x[ ] given on anN x N grid (nl, n2=0,l,...,N-l) by 
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Figure 5. (a) Quincunx lattice, (b) The corresponding Nyquist area in the frequency domain 
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If we now chain the down-sampling and up-sampling operators, we get 


W. 




\n 


\n 


_ Jxnj whenf^+n^isever 
Oelsewhere 


1 

2‘-" 


Xie^‘^) + Xie 


jiCO+TT) 


( 12 ) 


(13) 


Since quincunx sampling reduces the number of image samples by a factor of two, 
the corresponding reconstruction filterbank has two channels Figure 6. The low-pass filter H 
reduces the resolution by a factor of V2 ; the wavelet coefficients correspond to the output of the high-pass 
filter G [20-22], 


M M 



Figure 6. Perfect reconstruction filter bank on a quincunx lattice 


3.2.2. Fractional quincunx filters 

As starting point for our construction, we introduce a new 1 -D family of orthogonal filters; 




V 2 ( 2 -i- 2 cos<z>)^ 


(14) 


^(2 +2cos<z>)" +(2 —2cos<z>)" 

This is indexed by the continuously-varying order parameter a. Applying the diamond McClellan transform 

to the filter above is straightforward; it amounts to replacing cosco by 1 2 (1/ 2)(cos +COS (^2 ^ 

Thus, our quincunx refinement filter is given by [21-22]; 

a 

yf2(2 + COSCl\ +COSC02y 


HSen = 


^(2-|-C0SC9| -|-C0S®2)‘' +(2-C0SC9| - 00 ^, 002 )“ 


(15) 


This filter is guaranteed to be orthogonal because the McClellan transform has the property of 
preserving biorthogonality. Also, by construction, the a th order zero at co=7r gets mapped into 

a corresponding zero at (; this is precisely the condition that is required to get a 2-D wavelet 

H 

transform of order a . Also, note the isotropic behavior and the flatness of around the origin; i.e, 

) / V2=1 -h o(yr) V o 

II II Lxj ^ w jYiQ orthogonal wavelet filter is obtained by modulation 




( 16 ) 
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The corresponding orthogonal scaling function is defined implicitly as the solution of the quincunx two -scale 
relation; 


<p^ (x) = V2 hE i^x-n) 


(17) 


Since the refinement filter is orthogonal with respect to the quincunx lattice, it follows that 

)and that it is orthogonal to its integer translates. Moreover, for a >0, it will satisfy 
the partition of unity condition, which comes as a direct consequence of the vanishing of the filter at 


(^’^)Thus, we have the guarantee that our scheme will yield orthogonal wavelet bases of L 2 (R^). 
The underlying orthogonal quincunx wavelet is simply [23]; 


(x) = V2 22 hE ^Mx-n) 


(18) 


3.3. Spiht coding sheme 

SPIHT algorithm was propound by A. Said and Pearlman [24]. SPIHT performs a recursive 
partitioning to determine the position of the significant coefficients in the progeny of 
the considered coefficient. The basic principles of the SPIHT algorithm are a partial storage by amplitude of 
the wavelet coefficients, then a partitioning in hierarchical at each threshold applied, the trees are sorted on 
the basis of their meaning in two categories of trees and a scheduling of the transmission of the refinement bits. 
The following sets of coordinates are used to present the new coding method: 

O (i, j); set of coordinates of aU offspring of node (i, j). 

D (i, j); set of coordinates of aU descendants of the node (i, j). 

H; set of coordinates of all spatial orientation tree roots (nodes in the highest pyramid level). 

L (i, j)=D(i, j)-0(i, j); their algorithm stores the significant information in three ordered lists, 
to perform a nested coding, called insignificant sets list (LIS), insignificant pixels list (LIP), and significant 
pixels list (LSP), each entry is signified by (i, j) [23]. 

The coefficients insignificants in LIP are tested one more time to place those are significant in LSP, 
in the LIS, the coefficients that will be processed in the next refinement pass are found in the LSP list [25]. 
Eigure 7 shows parent-child relationship. 



Eigure 7. Parent-child relationship 


4. RESEARCH METHOD 

In image compression, minimize the size, keeping the interesting part and showing a good precision 
of the pathology remains an unresolved problem. This proposed method has been tested on several medical 
images to see its performance. According to LS method advantages, first of all, we focus on pathological 
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images, we segment the object of interest or the part that contains the disease, then after extracting the ROI, 
the different wavelets are applied on the new image without losing their quality of the image. 
Based on evaluation parameters such as PSNR and MSSIM, a comparison is given to determine the most 
powerful algorithm proposed. Figure 8 shows the proposed method steps. 



Figure 8. The proposed method steps 


4.1. Quality evaluation parameter 

4.1.1. The peak signal to noise ratio (PSNR) 

The peak signal to noise ratio (PSNR) used to measure the quality of reconstruction in ima^ 
compression. It refers to the ratio between signal and reconstruction error variance in Decibel scale. It can be 
represented as [26]: 


MSE = 


1 

MxN 


i = M 

I 


'fildJ) -/((,;)) 


(19) 


Where, MSB: mean squar error between two images. Mean Square Error (MSB) which requires two MxN 

grayscale images I and ^ where one of the images is considered as a compression of the other. 

The PSNR is defined as: 


PSNR= 10 log,. 


f 2 A 

{Dynamics of image) 

, MSE , 


( 20 ) 


4.1.2. The structural similarity index (SSIM) 

This parameter compares the similarity the brightness, contrast and structure between each pair of 
vectors, the structural similarity index between two signals x and y is given by the following 
expression [26, 27]: 

SSIM{x, y)=l{x, y)c{x, y)s{x, y) (2i) 


The MSSIM is used to measure the quality of the local image, which provides more information on 
the degradation of image quality, which is useful in medical imaging applications. It values exhibit greater 
consistency with the visual quality. 

1 M 

M55/M(/,/)=^£55/M(/,,/,) (22) 

^ i = l 

Where I and I are respectively the reference and degraded images, li and li are the contents of images at 
the i-th local window; M: the total number of local windows in image. Finally the quality measurement can 
provide a spatial map of the local image quality, which provides more information on the image quality 
degradation, which is useful in medical imaging applications. 
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5. RESULTS AND ANALYSIS 

The aim of our work lies in the possibility of reducing the rates, for which the image quality remains 
acceptable, which we are interested in hybrid compression methods based on activ contour segmentation and 
2D wavelet transforms because of their interesting properties, applied to medical images. 
Eor this reasons the PSNR evaluation parameters and the MSSIM similarity index are used to estimate and 
judge the quality of the compressed image. In order to verify the efficiency of our hybrid algorithm, 
we applied, (level set -i-CDE 9/7 (lifting scheme) -i-SPIHT) and (level set -i-quincunx -i-SPIHT) on MRI medical 
images of size 512x512 encoded by 8 bits per pixel. This image is taken from the GE Medical System 
database [28]. Eigure 9 illustrates the steps of the compressed image (region of interest image) 
by the quincunx wavelet. 



(a) (b) (c) 

Eigure 9. Image compressionsteps, (a) Brain original image, (b) The region of interest, (c) Quincunx 

wavelet decomposition 


The compressed image quality for different bit-rate values (number of bits per pixel) 
by the proposed algorithm is presented in the Eigure 10. Erom the results of the PSNR, MSSIM and MSE 
values obtained, we note that from 0.5bpp, the reconstructed image b ecomes almost perfect. 





Rc=0.25 PSNR=30.37 
MSSIM=0.7 


Rc=0.5 PSNR=38.09 
MSSIM=0.9 


Rc=0.75 PSNR=44.13 
MSSIM=0.97 





Rc=l PSNR=48.09 MSSIM=0.99 Rc=1.5 PSNR=57.22 MSSIM=1 


Rc=2 PSNR=66.32 MSSIM=1 


Eigure 10. Compressing of a ROI of an axial slice with level set coupled with quincunx wavelet and 

SPIHT coding 
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To ensure the performance of our proposed method, we compared between these different types of 
transform CDF9/7 (Lifting scheme)) coupled with level set method and SPIHT coding and CDF9/7 
(Lifting scheme)) coupled with SPIHT coding without applied the level set segmentation applied to the same 
medical image. The Table 1 and Table 2 given represent respectively these results. 


Table 1. PSNR andMSSIM variation using different methods proposed 


Type d’image 

RC(bpp) 

Level set + quincunx + SPIHT 

Level set + CDF9/7 (lifting) + SPIHT 

PSNR 

MSSIM 

MSB 

PSNR 

MSSIM 

MSB 


0.25 

30.37 

0.78 

55 

29.8 

0.75 

16.03 


0.5 

38.09 

0.9 

10.09 

36.14 

0.83 

7.01 

IRM 

0.75 

44.13 

0.97 

2.51 

42.06 

0.91 

2.5 

1 

48.99 

0.99 

0.82 

46.55 

0.96 

0.79 


1.5 

57.22 

1 

0.12 

57.26 

0,99 

0.11 


2 

66.32 

1 

0.02 

65.09 

0.99 

0.05 


Tab le 2. PSNR andMSSIM variation using different methods propo sed 


RC (bpp) 

Quincunx +SPIHT 

CDB9/7(lifting)+SPIHT 

PSNR 

MSSIM 

MSB 

PSNR 

MSSIM 

MSB 

0.25 

38.98 

0.91 

8.23 

36.35 

0.64 

15.22 

0.5 

44.33 

0.96 

2.4 

41.51 

0.79 

4.63 

0.75 

47.91 

0.98 

1.06 

45.88 

0.94 

1.7 

1 

50.75 

0.99 

0.56 

48.15 

0.96 

1.01 

1.5 

54.12 

0.99 

0.26 

53.04 

0.99 

0.33 

2 

56.45 

0.99 

0.16 

55.80 

0.99 

0.18 


Comparingthedifferent values of PSNR, MSSIM, andMSEwe ensuretheefficiency of our algorithm 
in terms of compressed image quality for the low bit-rate. The Figure 11 shows represent 
the results of the comparison obtained after applied the different algorithms proposed on an axial slice 
of brain imaging, given by the PSNR curve. 



-level set+Quincunx+SPIHT 

-level set+ CDF9/7(lifting scheme)+SPIHT 

-Quincunx+SPIHT 

-CDF9/7(lifting scheme)+ SPIHT 


1 1.2 
Rc(bpp) 



Figure 11. PSNR and MSSIM variation using different proposed methods 


The set of medical images in the following figures are used to ensure the performance of the proposed 
algorithm. For several images and from 0.5bpp our compression algorithm preserves the image structure and 
the pathologie form. It’s clearly shown the part, the type and the size of pathologie in all 
tested images. Figure 12 shows ultrasound imaging. Figure 13 shows sagittal section of brain. Figure 14 shows 
orthopedic imaging. 
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Original image 


PSNR=36.5, MSE=14.57 
MSSIM=0.93 
Quincunx-i-LS-i-SPIHT 


PSNR=33.12, MSE=31.7 
MSSIM=0.78 
CDP9/7-FLS-FSPIHT 


Eigure 12. Ultra sound imaging 



Original image 


PSNR=47.02, MSE=1.29, 
MSSIM=0.99 
Quincunx-i-LS-i-SPIHT 


PSNR=38.8, MSE=8.57, 
MSSIM=0.89 
CDE9/7-rLS-HSPIHT 


Eigure 13. Sagittal section of brain 



Original image 


PSNR=37.47, MSE=11.65, 
MSSIM=0.94 
Quincunx-i-LS-i-SPIHT 


PNR=33.65, MSE=30.63, 
MSSIM=0.79 
CDE9/7-rLS-rSPIHT 


Eigure 14. Orthopedic imaging 


According to the different results obtained for the several application of our hybrid method on color 
MRI medical image, we ensure that our proposed algorithm gives better quality and it is an efficient method 
compared with other techniques, more precisely when compressed the region of interest image after applying 
the LS segmentation with quincunx algorithm coupled with SPIHT coding, which gives a better performance 
and image quality-chapters. 
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6. CONCLUSION 

The work in this paper is to enhance medical images quality after the compression step. We used an 
hybrid method composed of the LS method and the quincunx wavelet transform applied on variety of medical 
images. We compared the founding results with the level set coupled with CDF9/7 lifting algorithm and SPIHT, 
we notice that the proposed algorithm gives better results than the other compression techniques. Generally 
the results obtained are very satisfactory in terms of compression ratio and compressed 
image quality. 
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